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ABSTRACT 

A wind passing over a surface may cause an instability in the surface such as the flapping seen when 
wind blows across a flag or waves when wind blows across water. We show that when a radially outflowing 
wind blows across a dense thin rotating disk, an initially flat disk is unstable to warping. When the 
wind is subsonic, the growth rate is dependent on the lift generated by the wind and the phase lag 
between the pressure perturbation and the vertical displacement in the disk caused by drag. When the 
wind is supersonic, the grow rate is primarily dependent on the form drag caused by the surface. While 
the radiative warping instability proposed by Pringle is promising for generating warps near luminous 
accreting objects, we expect the wind driven instability introduced here would dominate in objects which 
generate energetic outflows. 

Subject headings: instabilities; hydrodynamics; accretion disks 



1. INTRODUCTION 

For a wind blowing over a surface the velocity and pres- 
sure are related by Bernoulli's equation; P + \pv 2 is con- 
served along streamlines. Because the velocity increases as 
the wind passes over a protrusion of the surface, the pres- 
sure must decrease. This causes a force pulling the higher 
regions of the surface upwards. This force also causes lift 
on airplane wings. The Kelvin-Helmholtz instability oc- 
curs for the same physical reason at the boundary of two 
fluids of different densities moving respect to one another 
(e.g., Chandrasekhar 1961). A related instability exists 
for wind passing over water or for wind passing over fabric 
(e.g., Thwaites 1961; Mahon & Williams 1999). 

Near accreting compact objects such as active galac- 
tic nuclei, substantial amounts of kinetic energy can be 
present in a wind which may pass near denser colder ma- 
terial in the outer parts of an accretion disk. In this pa- 
per we consider the possibility that a wind passing over a 
dense disk can result in a warp instability similar to that 
caused by radiative forces (Pringle 1996; Maloney, Begel- 
man & Pringle 1996; Maloney, Begelman & Nowak 1997). 
Previous work, (Schandl & Meyer 1997; Porter 1998), has 
considered the torque on a disk which would be caused 
by ram pressure or the reaction force from a wind but has 
not explored the possibility that the wind/disk interaction 
could lift the disk, and so cause a warping instability. 



2. PERTURB ATIVE RESPONSE OF A RADIAL WIND TO A 
CORRUGATED SURFACE 

We divide the problem into two parts, a diffuse out- 
flowing wind and a dense infinitely thin disk. We first 
compute the effect of a vertical perturbation or ripple in 
the disk on the wind. The flow has a perturbation in the 
pressure along the surface of the disk which pushes on 
the disk. We then incorporate this force into the equa- 
tions of motion in the disk. This approach is similar to 
that used to investigate wind/ water wave interactions or 
wind/fabric interactions. We follow the perturbative ap- 



proach of Chandrasekhar (1961). 

We describe a warped surface by a displacement in the 
direction normal to its undisturbed plane, 



h(r, 6, t) = Re 



rjgi(rnO—k r r—uj€) 



(i) 



where S -C r. We assume a primarily radial flow with 
velocity u = uoe r + v where v is a perturbation, \v\ <C uq. 
We constrain the velocity so the component normal to the 
surface is zero; ^ — (Jj -I- u q-§^) h = v z . This constraint 
implies that on the surface 
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V • pv = 0, in cylindrical 



The continuity equation, 
coordinates is 

Id, . 19, d . dp . 

where p is the density in the wind. When the wind is sub- 
sonic, we search for a solution for the velocity and pressure 
perturbations in the wind that decays exponentially with 
increasing distance from the displaced surface. When the 
wind is supersonic, we search for a solution for the ve- 
locity and pressure perturbations that vary in phase with 
distance from the displaced surface. 

Pi = Re [pig{r, 0, z,i)] (4) 

v g =Re [v 6i ig(r,9,z,t)] 

iv = Re [v r ,ig(r,9,z,t)] 

v z =Re [v z> ig(r,9,z,t)] 

where 

g(r,6,Z,t) = e i{rnB-k T r-ut) 

and the Mach number, M = The pressure, P = 

Po+Pi, Pi <C Po and we adopt an equation of state P cx p T 
with sound speed in the wind c s . For the Kelvin-Helmholtz 
instability in an incompressible fluid k z is directly related 
to k x where the one fluid is moving with respect to the 



-kM for M < 1 
-ut+k,\z\) for m>1 
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other in the x direction. For incompressible potential flow 
solutions the vertical forms of the variables are solved ex- 
actly and they decay exponentially with height when the 
flow is subsonic and vary in phase with height when the 
flow is supersonic (e.g., Shivamoggi f998). 

In the the tight winding limit (k r 3> 1/V) and to first 
order the continuity equation becomes 



1 Pi / , 7 \ | i rnv e ,i 
— — (CJ + uok r ) + k r v r ,i — 



(6) 



= ik z v z ,isignz for M < 1 
= k z v z .isignz for M > 1. 
Eulcr's equation to first order in the same coordinate sys- 



tem 

V r ,t + U V ri r + UQ. r V r = 

u v e 

V0,t + U Ve, r H = 

r 

V z .t + U V Ztr = - 

Again in the tight winding limit 
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(8) 



(w + u k r )v z .i = ikzY^signz for M < 1 



= fc^f^signz for M > 1. 



From Equations (6) and (8) we find that 



(uj + u k r ) 2 +c 2 s (±k 2 z - k 2 



m 







(9) 



where the sign of the k 2 term is positive for M < 1 and 
negative for M > f. Neglecting the m 2 /r 2 term in the 
tight winding limit 

k 2 - + u k r ) 2 = k 2 for M < 1 (10) 

= -k 2 for M > f 

which relates the vertical exponential scale length or 
wavenumber to the radial wavenumber. 



When the wind is subsonic, the frequencies and wavevec- 
tors are real and k z > 0, then p\ < for z > 0. Where 
the surface is high, above the surface the pressure is lower 
than the mean and below the surface the pressure is higher. 
This causes a destabilizing force which we also call lift in 
analogy to the aerodynamics of wings. Because the sign of 
the pressure perturbation is opposite on either size of the 
disk (oc sign z), the pressure differential across the sur- 
face exerts a force on the surface with magnitude twice 
Pi in the direction normal to the surface (see Figure f). 
The pressure differential across the surface caused by the 
wind should exert a force on the surface. An instability 
should exit when m — which is a direct analogy for the 
Kelvin-Helmholtz instability or for the instability of a wind 
passing over fabric. 

When the wind is supersonic, the pressure perturbations 
are 90 degrees out of phase with the surface perturbations. 
The phase of the pressure perturbations should increase 
with both r and | z | which implies that the sign of k z is the 
same as the sign of k r . Because the sign of the pressure 
perturbation is opposite on either size of the disk, there 
is a drag force on the surface which is not present when 
the flow is subsonic and approximated as a potential flow 
solution. Because shocks should form we expect the actual 
flow to be more complicated than given by the above equa- 
tions (see Figure 2 for an illustration). When the shocks 
and expansion waves lag the high points of the surface, we 
expect a vertical force on the high points of the surface 
which pushes the surface towards the midplanc. When 
m = 0, instead of a destabilizing lift force, we expect a 
stabilizing vertical force. But since there is drag, energy 
from the wind is transferred to the surface and the ampli- 
tude of the surface perturbation should increase (e.g., §5.5 
of Shivamoggi f998 on the Kelvin Helmholtz instability 
with a supersonic flow). 
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Fig. 1— As a subsonic wind passes over corrcgations, 
pressure changes cause lift at the high points of the 
surface. 

Using Equations (2, 8) we can see that 

|i = -§£(w + u k r ) 2 signz for M < 1 (11) 
= -^JL(tj + u a k r ) 2 forM>l 



Fig. 2— When the wind is supersonic, the maxima 
of the pressure perturbations vary as a function of 
distance from the surface. Shocks cause pressure in- 
crease, and expansion fans cause pressure decreases. 
The pressure differences across the surface result in a 
drag on the surface. When the expansion waves and 
shocks lag the high points of the surface, there is a 
vertical force on the surface which pushes the high 
points towards the midplanc. 

When m = 1 for both subsonic and supersonic flows we 
must consider the torque on an annulus, which is the cross 
product of the radial vector and the pressure differential, 
integrated azimuthally about the ring. 

3. THE TORQUE ON DISK ANNULI 

To integrate the torque azimuthally we now shift nota- 
tion and follow that used by Pringle (1996) to describe the 
warped disk. The tilt vector for a ring at radius r 

I = (cos 7 sin 0, sin 7 sin 0, cos 0) (12) 
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where (3(r, t) is the local angle of tilt of the disk with re- 
spect to the z axis, and the descending node of the disk 
material is at an angle j(r, i) — 7r/2 to the x axis. The unit 
vector toward points on the surface 

i(r, 4>) — (cos <p sin 7 + sin <p cos 7 cos (3, (13) 
sin 4> sin 7 cos (3 — cos <p cos 7, — sin <fi sin (3) 

where (f> = at the descending node and 7 and /3 are 
both functions of r. The external coordinate system, 
x = (x,y,z), which in cylindrical coordinates is (r, 9, z) 
where 9 = tan _1 (j//x), can be related to that described 
by Equation (13). When (3 < 1 

6» = ^ + 7-tt/2 (14) 

and the vertical displacement of the surface or z compo- 
nent of x 

h(r, 9, t) = ~(3r sin <j> = -f3r cos(6> - 7) (15) 

To relate this formalism to that used in the previous sec- 
tion we set 7 = k r r + ut, m = 1 and S = —fir so that 

h(r,6,t) = -f3rcos(6-k r r-ujt) (16) 
which is in the same form as Equation(l). 

Pringle (1996) defines an elemental area vector dS — 
sT- x s^drd(j). where s" r = |^ and = The normal to 

the disk surface n = . To first order in (3 

\dS\ 

s ^ = / x x = (cos (4> + 7) , sin (<fi + 7) , — /3 cos 0) (17) 

h = I — x(j3r-f' cos (f) ~ r (3' sin (/)) (18) 

(see Pringle 1996, Equations 2.11 and 2.17) where /?' and 
7' refer to the derivatives of (3 and 7 with respect to r. 
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Fig. 3- The lift on an annulus causes the annulus to 
precess. When the flow is supersonic, the force on the 
high points of the surface is towards the plane and the 
annulus should precess in the opposite direction. 

When the wind is subsonic, in the previous section we 
found to first order that the pressure differential across the 
surface was in phase with the corrugations of the surface. 
However, the pressure in the region of laminar flow should 
actually be somewhat lower on the leeward side of each 
trough and so that there would be a lag between the pres- 
sure and the surface variations (Jeffreys 1924). Kendall 
(1970) measured sinusoidal pressure variations in response 
to a wind passing over a sinusoidally varying ruber surface 
and showed that the pressure was indeed offset in phase 
with the surface. We assume that the pressure difference 
across the disk surface can be described 

AP = 2 Pl sin(0 + S) forM<l (19) 



where the phase lag between the pressure and the surface 
is given by <5. For k r > we expect S > and that S < 
for k r < 0. 

When the wind is supersonic, shocks should cause an 
effective phase shift in the pressure perturbation on the 
surface. In this case 

AP = 2p[ cos((/> + S) for M > 1 (20) 

where 

Pi = 5£ ^(u 2 + u k r f (21) 

see Equation(ll). The sign of p[ depends on the sign of 
k ri since the sign of k z depends on the sign of k r . 

The pressure differential across the surface exerts a force 
on the disk parallel to the normal to the surface, fi. The 
resulting torque per unit mass per unit mass for an annulus 
is the integral 

1 f 27r AP 
t w = ^ J rax n—d<f> (22) 

where X is the disk surface density (mass per unit area). 
Using the previous 4 equations, we perform this integral 
for small (3 and S 

tZ, — ^-[( sm 7 ; — C0S7, 0) - <5(cos7, SU17, 0) (23) 
+ (0,0,(5/3)] forM<l 

= ?jf{— <5(sin7, — cos 7, 0) — (cos 7, sin 7, 0) 
+(0,0,/?)] forM>l 

To first order in (3 the angular momentum axis of an an- 
nulus is I = (/3 cos 7, /3sin7, 1). From comparing Equation 

(23) to % we can see that the first term in this equation 
causes the disk to precess, and the second term causes the 
tilt to increase or decrease. The tilt increases when 5 < 
or k r < 0, for a subsonic wind and when p[ < or k r < 
for a supersonic wind (see Figures 3 and 4). The third 
term corresponds to a rate of change in the total angu- 
lar momentum of the annulus which would cause a small 
amount of outflow in the disk. 




Fig. 4- The drag on the surface causes a vertical force 
on the disk. This results in a torque which can either 
increase or decrease the tilt. A ring is shown tilted 
towards the line of sight. The slope of the warp with 
respect to the wind depends upon whether the warp is 
leading or trailing. Consequently the torque either in- 
creases or decreases the tilt depending upon whether 
the warp is leading or trailing. 
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3.1. Equation of motion for the annulus 
The angular momentum per unit mass of an annulus of 



radius 



~VL where fl is the angular rotation rate. 



When radial motions in the disk are small 

o „ dl , 
r 2 n— = Tg + T v + T w (24) 

where f g is the torque from gravity (equal to zero when 
the potential is spherical) and f v is the torque due to vis- 
cous forces. As done by Pringle (1996), we set W = fie 11 
so that I = (Re[W], Im(W], 1). Ignoring the z component 
of the torque and using Equation(23) we find 



w 



for M < 1 
for M > 1 



(25) 



where we have ignored gravitation and viscous terms. Us- 
ing Equations(10, 11, 21) for pi and pi, 

dW rp \k r \M 2 f . x ^ w w M ^ 1 (2g) 



dt 



rP°*rM 2 _ ^ + w for M > 1 



where we have assumed slowly varying modes, w <C uo^r- 
We define a parameter Q which describes the ratio of the 
kinetic energy in the wind to the rotational energy in the 

disk, Q = ~ P ^"° where p w is the density of the 

wind and v c = fir is the velocity of a particle in a circular 
orbit. We can write 

dW^_ gmr r B(i + s)w forM<1 (2?) 



at 



Vl-M 2 

+ forM>l 



(28) 



The precession rate 

^ w w « -Q|fc r r|ft for M -C 1 

«^Q|fc r r|0/Mfor M> 1 

where to simplify the expression we have assumed M -C 1 
or M > 1 and ignored the dependence on the adiabatic 
index. 

For a subsonic wind, growing modes occur when S < 0. 
We expect that the phase lag S < when fc r < so that 
growing modes will occur for 7' < 0. The growing mode 
will have lines of nodes following a trailing spiral instead of 
a leading spiral as was the case for the radiatively driven 
warp instability (Pringle 1996). For a supersonic wind, 
growing modes occur when k r > and so will have a lead- 
ing spiral shape. 

We expect the growing mode to grow as (3 oc e*/ tm with 
growth timescale 

t w ~± forM<l (29) 
- £ for M > 1. 



n - 



In terms of the dynamical time, td 

u 

for M > 1 



Q^rs fOT M < 1 

M 



(30) 



When the wind is supersonic, the growth rate exceeds the 
precession rate and the amplitude of the perturbation can 



grow quickly. We see that the precession and growth rates 
depend on Q, the ratio of the kinetic energy in the wind 
and the rotational energy in the disk. Though we have 
neglected the viscous forces (in Equation 22), we expect 
them to damp the growth of unstable modes primarily 
for the short wavelength modes. We expect the disk self 
gravity to affect the precession frequency but not strongly 
affect the growth rate of the mode. 

Since we have followed the notation by Pringle (1996), 
we can directly compare the importance of the torque 
caused by the wind to that caused by the absorption 
of radiation from a central source. Ignoring the viscous 
damping term, the precession rate caused by the radiative 
torque is given by uj t = 12 J^^n c (Equations 3.6, 3.8 by 
Pringle 1996). Comparing this to our Equations(27,28) we 
find that the wind driven warping instability is likely to 
dominate when 

Jj-<P M*r*~p w uy. (31) 

We can interpret this inequality in terms of the momentum 
flux. When the momentum flux from the wind dominates 
that from a central radiation source, a wind driven insta- 
bility may dominate. This suggests that accreting sources 
which impart more energy into driving winds than into 
radiation would be more likely to drive wind driven insta- 
bilities in their accretion disks. 

4. THE PHASE LAG AND THE CRITICAL ANGLE FOR 
BOUNDARY LAYER SEPARATION 

The simplest description of subsonic wind flow near a 
corrugated surface such as a flag or water wave is an irrota- 
tional flow above the interface with wind velocity averaged 
over one wavelength that is constant with height. This is 
a "potential flow" (Lamb 1932) and since the pressure is 
related to the wind velocity via Bernoulli's equation, it is 
in exact anti-phase with the surface (as derived perturba- 
tively above). So the energy flux between the wind and 
wave is zero. For the wind to do work on the wave, the 
pressure must be shifted in phase relative to the potential 
solution. Jeffreys (1924, 1925) suggested that the down- 
wind side of the wave is sheltered from the wind so that 
the pressure on this slope is reduced there but increased 
on the upwind side. This causes a drag force, and as a 
result waves traveling with the wind increase in amplitude 
and those traveling against the wind are damped. Scal- 
ing based on idealized theory for this (Miles 1957, and 
subsequent work) unfortunately does not predict the ex- 
pcrimentaly measured growth rates very well (e.g. Kendall 
1970; Riley, Donelan & Hui 1982; Donclan 1999). 

In the theory outlined by Miles (1957), the phase lag 
is set by numerical constants determined from turbulent 
mean profiles and does not depend upon the amplitude of 
the surface variations. We can understand this consider- 
ing the difference in the rate that a turbulent boundary 
layer grows on the leeward and windward sides of a wavy 
surface. Because the growth rate of the boundary layer 
depends on the pressure gradient along the direction of 
flow outside the boundary layer, we expect a difference in 
the boundary layer thickness on the leeward and windward 
sides that is proportional to the pressure difference across 
the wave, and this in turn is proportional to the amplitude 
of the wave. However, the thickness of the boundary layer 
may be set by a turbulent eddy viscosity and so should 
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also be proportional to the amplitude of the surface per- 
turbation. Thus, we predict that the phase angle is not be 
strongly dependent on the amplitude of the surface varia- 
tions. 

In most astrophysical situations we expect high 
Reynolds number flows and so can draw on the theory 
of aerodynamics. A laminar or viscous boundary layer 
next to a surface is only stable when the pressure gradi- 
ent along the surface in the direction of flow is negative. 
However this condition is violated to the leeward of a cor- 
rugation. The same situation occurs along airplane wings 
and the boundary layer becomes turbulent. Because of the 
Coanda effect, turbulent boundary layers arc less likely to 
separate from the surface than viscous boundary layers. 
This allows the boundary layer to remain near the surface 
along airplane wings even though the pressure gradient 
in the laminar flow region is not favorable. For most air- 
plane wings the boundary layer does not separate from the 
surface until the angle of attack from the wing is > 15°, 
though this critical angle can increase to 45° for short thick 
wings. Past this critical angle, little or no lift is generated 
and the wing 'stalls'. We expect a similar effect in our 
astrophysical disk. When the amplitude of the corruga- 
tions reaches a critical slope, (a critical value for \Sk r \), 
we expect that the boundary layer will separate from the 
surface, the lift will decrease causing the precession rate to 
drop, and the drag will increase. At this point we expect 
the mode to saturate. 

When the wind is supersonic, a smooth continuous flow 
is unlikely (see Figure 2). A full calculation would re- 
quire resolving shocks and expansion flows. However the 
perturbative method we have used should allow us to at 
least estimate the magnitude of the instability. We do ex- 
pect a major difference in the character of the wind flow 
near the surface when the slope of the surface exceeds the 
Mach angle, or when \Sk r \ > atanM^ 1 . Past this angle 
a turbulent boundary layer should develop and the insta- 
bility may saturate. In this regime it may be possible to 
estimate the torque on the disk using an approximation 
(that of Newtonian flow) which primarily considers ram 
pressure or reactive force on the disk surface (e.g., Porter 
1998; Schandl & Meyer 1997). 

5. DISCUSSION AND SUMMARY 

In this paper we have outlined a possible instability 
which could occur in accreting objects with energetic out- 
flows. A wind passing over a rotating thin dense disk is 
likely to cause a warping instability to grow in the disk. 
This instability could occur in situations where an ener- 
getic outflow driven by a compact object passes over a 
dense disk, such as might happen in binary X-ray sources 
or active galactic nuclei. In active galactic nuclei, at radii 
of order a parsec from a black hole there is evidence for 
the existence of dense warped disks from maser observa- 
tions (e.g., Herrnstein et al. 1996). Large scale outflows 
are predicted for a variety of types of accretion flows and 
seem to be an integral part of these flows. Observational 
evidence for disk winds is fairly ubiquitous and seen in 
both X-ray binaries (Brandt & Schulz 2000; Chiang 2001) 
and in AGNs (e.g., Murray & Chiang 1995). 



When a subsonic wind passes over a corrugated dense 
disk, the corrugations will cause pressure variations along 
the surface of the disk. These pressure variations cause 
lift, which cause the annuli in the disk to precess. Form 
drag caused by the corrugated surface will cause pressure 
variations in the disk to lag the vertical displacement of the 
disk. This causes a torque on the warped disk which can 
increase the amplitude of the perturbation. We expect this 
instability to saturate or cease growing at a critical slope 
when the lift decreases due to boundary layer separation. 

When a supersonic wind passes over a corrugated disk, 
the disk primarily experiences form drag. The drag causes 
a torque on the warped disk which can increase the am- 
plitude of the corrugations. When there is an effective 
lag between the pressure and the vertical velocity of the 
surface, there will be a vertical force on the disk which 
causes it to precess. We expect a major change in the na- 
ture of the flow when the slope of the corrugations exceeds 
the Mach angle. While we primarily expect accreting as- 
trophysical objects to drive supersonic outflows, if a thick 
outflowing subsonic boundary layer is created, the sub- 
sonic theory outlined here may be appropriate inside this 
layer. The instability can dominate the radiative induced 
warping instability in sources where the energy released in 
an outflow (either sub or supersonic) exceeds that emitted 
as radiation. 

In this paper we have concentrated on a purely hydro- 
dynamic flow. However we expect the magnetic field to be 
dynamically important both in the disk and wind. If the 
wind and disk are magnetized, flux freezing would make 
the wind and disk act as mutually impenetrable bodies. 
The wind would be less likely to cause hydrodynamical 
instabilities that could destroy the disk. In this paper we 
have focused on the affect of a radially outflowing wind on 
the planarity of a disk in the absence of radiation. Future 
investigations could develop a prediction for the phase lag 
angle for subsonic wind flows (a historically difficult hy- 
drodynamics problem), and calculate the details of the 
shocks which are likely to be present when the wind is su- 
personic. In this paper we have restricted our study to 
radially outflowing winds, however astrophysical outflows 
typically have a nonzero, radius dependent, azimuthal ve- 
locity. We have also neglected the role of gravitational, 
magnetic and viscous forces. Future work can consider 
the role of these forces, include rotating winds and explore 
the effects of both radiation and winds on the planarity 
of the disk, particularly in the regime where the winds 
themselves are radiatively driven. 
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sities for Research in Astronomy, Incorporated, under 
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